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ABSTRACT 

We have conducted interferometric observations with the Combined Array for Research in Milhme- 
ter Astronomy (CARMA) and an on-the-fly mapping with the 45-m telescope at Nobeyama Radio 
Observatory (NR045) in the CO {J— 1-0) emission hne of the nearby spiral galaxy NGC 3521. Using 
the new combined CARMA -I- NR045 data of NGC 3521, together with similar data for NGC 5194 
(M51a) and archival SINGS Ha, 24 ^m, THINGS H i and GALEX FUV data for these two galax- 
ies, we investigate the empirical scaling law that connects the surface density of star formation rate 
(SFR) and cold gas (known as the Schmidt-Kennicutt law or S-K law) on a spatially-resolved basis, 
and find a super-linear slope for the S-K law when carefully subtracting the background emissions 
in the SFR image. We argue that plausibly deriving SFR maps of nearby galaxies requires the dif- 
fuse stellar and dust background emission to be subtracted carefully (especially in the mid-infrared 
and to a lesser extent in the FUV). Applying this approach, we perform a pixel-by-pixel analysis on 
both galaxies and quantitatively show that the controversial results whether the molecular S-K law 
(expressed as Ssfr oc S^^^ ) is super-linear or basically linear is a result of removing or preserving 
the local background. In both galaxies, the power index of the molecular S-K law is super-linear 
(7H2 ^ 1-5) at the highest available resolution (^^230 pc), and decreases monotonically for decreasing 
resolution. We also find in both galaxies that the scatter of the molecular S-K law ((TH2) monoton- 
ically increases as the resolution becomes higher, indicating a trend for which the S-K law breaks 
down below some scale. Both 7H2 and are systematically larger in M51a than in NGC 3521, but 
when plotted against the de-projected scale ((5dp), both quantities become highly consistent for the 
two galaxies, tentatively suggesting that the sub-kpc molecular S-K law in spiral galaxies depends 
only on the scale being considered, without varying amongst spiral galaxies. A logarithmic function 
7H2 = —1-1 log[(5dp/kpc] + 1.4 and a linear relation tTH2 = —0.2 [(5dp/kpc] + 0.7 are obtained through 
fitting to the M51a data, which describes both galaxies impressively well on sub-kpc scales. A larger 
sample of galaxies with better sensitivity, resolution and broader field of view are required to test the 
general applicability of these relations. 

Subject headings: galaxies: individual (NGC 3521, NGC 5194, M51a); galaxies: ISM; galaxies: spiral; 
ISM: molecules; radio lines: galaxies. Stars: Formation 



1. INTRODUCTION 

Much of the information we obtain on galaxies and 
their evolution is the consequence of the formation of 
stars, the process that depletes galaxies of their gas, dy- 
namically energizes the stellar disks and enriches the in- 
terstellar medium with heavy elements through genera- 
tions of evolving stars. Despite its fundamental impor- 
tance on nearly all scales of astronomy, star formation is 
an exceedingly complex process to explore in the ab initio 
manner as a result of the many entangling aspects includ- 
ing the dynamical, geometrical, magnetic and chemical 
properties of ISM, as well as the feedback that the young 
stellar objects exert on their surrounding environment by 
radiation, stellar wind and explosion of supernovae. 

In spite of this difficulty, the well-established empirical 
scaling laws of star formation that connect the large- 
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scale star formation rate (SFR) to interstellar medium 
(ISM) properties already play a crucial role in establish- 
ing plausible scenarios for galaxy evolution, and serve as 
essential prescription in modeling and simulating galaxy 
formation and evolution. Among them, the most widely 
adopted power -law relationship known as Schmidt law 
(jSchmidtl 119590 links t he volume or sur face density of 
SFR to that of the gas (jKennicutt 1 11989^). From an ob- 
servational point of view, a variety of SFR indicators are 
currently used, along with gas density information de- 
rived from emission lines of molecules (e.g. CO, HCN), 
to trace the relation between the birthrate of stars and 
their environment. By averaging over the whole disk of 
individual galaxies and taking the total (atomic -I- molec- 
ular) cold gas into account, a single power-law scaling 
between the surface density of star formation and gas 
(the Schmidt-Ke nnicutt law, abbreviated as the S-K law; 
iKennicutt 111 9981 ) extending over five orders of magnitude 
in the gas density has been well established, taking the 
form SsFR oc Ei,4±0-i5 where EgFR (in Mq yr"! kpc-^) 
and Sgas (in M© pc~^) denote the disk-averaged SFR 
and gas surface density, respectively. 

Both parameters in this equation are, however, global 
quantities averaging over local values that may vary 
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by several orders of magnitude. Further understand- 
ing of star formation on a more physical basis requires 
scrutinizing galaxies at the sub-galactic scale. Lim- 
ited by development of the relevant instrumentation, 
spatially-resolved investigations of the gas content of 
individual star-forming sites have become a front-line 
topic only in recent years (e.g . . IWoiig fc B litz 2002; 
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2007; Bigicl et al. 2008!: iBlanc et al.l 120091 ). Spatiallv- 
resolved studies promise to yield insights into the physi- 
cal processes driving the scaling laws, much more directly 
than the disk-averaged results. Measuring the star for- 
mation and gas at sub-kpc scales requires high resolu- 
tion and sensitivity, which are currently available only 
for small or limited sample s both for ato mic hydrogen 
(e.g., through THINGS; Walter et al.l l2008D or molecula r 
hydrogen (via CO, e.g., HERACLES: iLerov et al.l [2009h. 
Observational studies on the slope of the local S-K law in 
a handful of nearby galax ies result in power-law indices 
of 1- 3 (reviewed by, e.g.. lElmegreen I [20021 : iBigiel et al.l 
120081) ■ where the cold gas mass in both atomic and molec- 
ular phases are taken into account. However, in recent 
times, it has become clearer that star formation activ- 
ity, at least on local scales, is much less co r related with 
Hi than with CO ([Kennicutt et al.l l2007t IBigiel et al.l 
120081 hereafter K07 and BOS, respectively), which sug- 
gests that further investigations on the power law index 
of the S-K law should concentrate on the scaling rela- 
tionship between SFR and molecular gas. 

The availability of high resolution and sensitivity maps 
for CO as the most commonly used molecular hydrogen 
tracer thus provide the main actual limitation to sub- 
kpc studies of star formation scaling laws. In general, 
interferometric CO maps have higher resolution (typi- 
cally a few arcseconds) and are powerful in resolving fine 
structures, but their limited field of view (FoV) and the 
shortcoming of missing the zero spacing fiux often down- 
grade their scientific value; on the other hand, despite 
that single-dish observations are time-saving (especially 
if On-The-Fly mapping techniques are introduced) and 
cover much wider FoVs, their larger beams often smear 
out the details of star forming sites. In recent years, along 
with use of mosaicing that improves the FoV of interfer- 
ometry, the technique of compensating the absent short 
baselines components with single-dish data (at the price 
of losing some resolution) has been developed. With this 
technique mature, single-dish-compensated interferomet- 
ric data of CO transition line maps become the best 
available molecular gas tracer for investigations of star 
formation scaling laws. In this direction, a new CARMA 
(the Combined Array for Research in Millimeter-wave 
Astronomy) & NR045 (Nobeyama Radio Observatory 
45m telescope) CO surv ey of nearby galaxies is under- 
way (jKoda et al.ll2009bl ) with the goal of better under- 
standing the evolution of the ISM, galactic dynamics, 
and the star formation laws. Included in the sample of 
this CO survey, one of our program galaxies, NGC 3521, 
is mapped by CARMA with a 3.65"x3.14" synthesized 
beam and by NR045 using a 18.4" beam, leading to a 
final resolution of 5.43" x 4.68" when the two sets of data 
are comb ined in uv-space wit h the new technique estab- 
lished bv lKodaet all ([2009al) . This CARMA -I- NR045 
survey not only isolates typical GMCs with its ^220 pc 
resolution at a distance of ~8 Mpc (a typical CMC sep- 



aration is ~200 pc in the Milky Way, iKoda et al.ll2006D . 
but also detects the interarm extended emission at a high 
sensitivity (Icr ^ 3 Mq pc~^). 

The other crucial parameter in the resolved S-K law 
study is the SFR. Over the past twenty to thirty years, 
extensive effor ts have b e en rn ade to plausibly derive 
SFRs (see Calzetti et all (l2qIol) for a review of the last 
decade and lKennicutt I ( 19981 ) for earlier times). Conven- 
tional optical SFR tracers like Ha often severely suffer 
from dust extinctio n (Ay ^ 2.2 mag i n typical extra- 
galactic H II regions. fCalzetti et al.l[2007l ). which changes 
drastically from location to location and can be com- 
pletely obscured in dense star forming re gions where Ay 
can reach ^6 mag (|Scoville et al.ll200l for the case of 
M51a). Recent space-borne facilities have opened in- 
frared (IR) and ultraviolet (UV) windows for star forma- 
tion stu dy. The Spitzer I nfrared Nearby Galaxies Survey 
(SINGS, EennicutLet al. 200 3) and the GALE X Nearby 
Galaxies Survey (NGS, Bia nchi et al.l l2003al lH) enables 
us to image the details of dust obscured and unobscured 
star formation sites in parallel, so that both the dust 
reprocessed and directly emerging radiation from young 
stellar associations are measured. K07 and Calzetti et al] 
(j2007| ) have justified the feasibility of correcting the num- 
ber of ionizing photons (as traced by the Ha recombina- 
tion line) for the effects of dust extinction by adding a 
weighted component from the Spitzer MIPS 24 fj,m lu- 
minosity in individu al star formi ng regions. This Ha-|-24 
/xm SFR motivated ILerov et aH (2008) and B08 to pro- 
pose another composite SFR tracer that corrects dust at- 
tenuation of far-UV (FUV) surface brightness using the 
same mid-IR component (but with different weights). 

The extinction-corrected Ha line emission has long 
been proved a reliable estimator of the rate of forma- 
tion of young, massive stars, since stellar synthesis mod- 
els show that for co-eval stellar populations with stan- 
dard Initial Mass Functions that rate remains fairly con- 
stant for the first few million years (~5 Myr), and then 
abru ptly fades away as massive stars become supernovae 
f Scha erer fc de Kot er 1997). 

In marked contrast, a stellar population can emit in 
the continuum UV (longward of 1300 A) for up to one 
order of magnitude longer timescale. The Local Vol- 
ume Legacy (LVL) Survey team has found many dwarf 
galaxies that show little or no sign of dust and sporadic 
highly concentrated Ha emission, yet show strong and 
extended UV emission and have red UV colors, consis- 
tent with st ellar p opulations at the age of a few 100 Myrs 
pale et al.l 120091) . This comphcates the comparison be- 
tween Ha-based and UV-based SFR indicators. An in- 
depth discussion on the consistency and discrepa ncy of 
derivi ng SFR via Ha and UV is presented in Le e et al.l 
((20091) . 

The issue of a molecular-only S-K law is only now start- 
ing to be tackled. Several earlier attempts were made 
by utilizing azimuthally averaged data. Employing this 
strategy, Wong fc Bhtz (2002) studied seven CO-bright 
spiral galaxies and concluded that SFR is in dire c t pro - 
portion to molecular cloud density; Boissier et al.l ()2003f ) 
compiled the data for 16 galaxies and found the slope 
~0.6-1.3 (though most data points seemi ngly favor a 
super- linear slope; see Figure 10 therein); iHever et al.l 
(j2004| ) examined the low-luminosity, molecule-poor M33 
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galaxy by combining 60 and 100 fim IR images with CO 
(J=l-0) single-dish line map, and suggested a power in- 
dex of 1.36±0.08. More recently, SFR and molecular gas 
surface density have been compared on individual sub- 
kpc regions. K07 performed an analysis of star formation 
complexes in M51a (NGC 5194) and obtained a slope of 
1.37±0.03 for the molecular S-K law at ~500 pc reso- 
lution. On the other hand, BOS performed a pixel-by- 
pixel analysis (750 pc scale) for seven spiral galaxies and 
found a power index 1.0±0.2 or a linear molecular S-K 
law. Thus it can be seen that the results from different 
groups as yet fall into two categories: linear or super- 
linear scaling r elatio n, which is still now under debate. 
iWong fc Blitd (I2002D wo uld have obtained results simi- 
lar to lHever et al.l ()2004[ ) if radially dependent extinction 
corrections were applied to their Ha data. BOS qualita- 
tively suggests that the inconsistency between their re- 
sults and K07 is likely an effect of removal of the local 
background emission in the K07 photometry of the SFR 
tracers. Clearly, background subtractions in the SFR 
maps is a crucial component for the investigation of star 
formation laws. 

In parallel to these observational efforts, various star 
formation models have been constructed, generally as- 
suming one mechanism that dominates the formation 
process. Among these models, a postulated constant 
star formation efficiency independe nt of the environ - 
ment will lead to a linear relation (jLerov et al.l 120081 ) . 
in which case the mere presence of the gas is suffi- 
cient for star formation; a slope of 1.5 will be found 
if star formation is driven by large scale gravitational 
instabilities, and thus occurs in a dynamical ( free-fall) 
timcscale of molecular gas (|Elmegreen I l2002f ): cloud- 
cloud collisions will drive the slope to be as steep as 
~2 (jTask cr fc Tan 2009i): other cons idered mechanisms 
includ e global resonances (|Wvse fc S ilk IQS j, galactic 
shear (|Hunter et al.lll998D. dust s hielding (Gnedin et al.l 
I2009D. turbu lence (iKrumholz et a l. 2009) and others (for 
reviews, s ee lElmegreen I l2002t iMcKee fc Ostrikeii 120071: 
E3[2Qp). 

In this paper, we present a study of the sub-kpc S-K 
law in two nearby galaxies, NGC 3521 and M51a (NGC 
5194). K07 and B08, despite sharing data sets for M51a, 
obtain different power law indices (super-linear vs. lin- 
ear) of the molecular-only and total-hydrogen S-K laws. 
By employing a new method of local background subtrac- 
tion, we show quantitatively that including or removing 
the local background in the maps tracing the SFR has a 
profound effect on the final results in M51a, as discussed 
in BOS. By analyzing our high quality data, we infer that 
the linear correlation found by BOS for the molecular S- 
K law is a consequence of inclusion of stellar local back- 
ground unrelated to current star formation. Throughout 
this paper, we base our analysis on the Ha-|-24 ^m ap- 
proach of calibrating SFRs, albeit the FUV-h24/^m SFR 
is included for the purpose of comparison. 

Among the extensive measurements of the distance to 
these two nearby galaxies, we adopt the values from the 
compilation by the LVL survey group: S.03 Mpc for NGC 
3521 (from the flow- field correct ed recessional velocity , 
assuming Ho=75 km s"^ Mpc~^: [Kennicutt et"aI1l200Sf ). 
and S.OO Mpc for M51a (by direct distance m easurement 
of the M51 group, iKarachentsev et al.l[200l . 



NGC 3521 and M51a possess similar properties in sev- 
eral aspects. First of all, the disk-average SFR sur- 
face density is -0.0026 Mq yr'^ kpc'^ and ~0.0057 
M(n yr"-'^ ^P*?."^ these two spirals, respectively 
(Calzetti et al.l[2bl0t ). characterizing both as quiescently 
star-forming systems. Moreover, both galaxies are 
shown to be metal-rich from spectroscopic measure- 
ments, with 12-Hog (0/H) values 9.01 for NGC 3521 
and 9.18 for M51a (|Moustakas et al.l l2010l ). In addi- 
tion, cl assified as SABbc in the Hubble sequence (RC3, 
Ide Vau coulcurs et al. 1 99lb. NGC 3521, with an inclina- 
tion of 72.7° ( de Blok et al.l 120081 ). is a flocculent-type 
spiral g alaxy that has a tightly wound two-arm pattern 
(|Thorn lev 199( |), while its muc h less inclined peer M51a 
(z = 42°, Tam burro et al.ir200Sf ). a grand-design spiral of 
morphological type SAbc, shows a prominent two-arm 
configuration. The virtually identical size of their galac- 
tic disks is a consequence of having almost identical dis- 
tances (~S Mpc) and apparent size in the Ha images 
that depicts unambiguously the spiral arms (284" and 
296" in radius, respectively). These facts, when summed 
together, suggest that the two star-forming disks are 
very similar, with the major difference being their incli- 
nations. This is an important ingredient for understand- 
ing the influence of spatial resolution on the resultant 
shape of star formation laws (see the result and discus- 
sion section). 

The paper is organized as follows: section 2 describes 
the data and, especially for the CARMA/NR045 data, 
provides some detail on the observations and data reduc- 
tion; section 3 presents our new method of automated lo- 
cal background removal which enables our pixel-by-pixel 
analysis to reproduce both the results of K07 and BOS ac- 
curately, and discusses the necessity of local background 
removal; section 4 is devoted to applying our strategy 
to both spirals to test the behavior of the sub-kpc star 
formation laws; the results are discussed in section 5 and 
summarized in section 6. 

2. DATA 

In Figure [1] we show the Ha (KPNO 2.1-m telescope), 
24 fim (Spitzer SINGS survey), CO (our CARMA ob- 
servation) and FUV (GALEX NGS survey) images for 
NGC 3521. More technical details of these data and the 
CARMA -I- NR045 CO map of M51a are given in the 
following subsections. The Ha, 24 /im, BIMA SONG 
CO and FUV images of M51a, shared with K07 and BOS, 
have been thouroughly described in those papers and de- 
tails will only be provided when necessary. To compare 
our analysis with those in previous studies, we use the 
BIMA SONG data which were used in those studies, but 
our CARMA 4- NR045 data are employed for further 
investigations. 

2.1. CARMA and NR045 CO (.1=1-0) Data 

The CO ( J=l-0) interferometric observations of NGC 
3521 were undertaken using the Combined Array for Re- 
search in Millimeter Astronomy (CARMA) from Febru- 
ary to March in 2009. CARMA is the combination of 
the six 10 m antennas of the Owens Valley Radio Obser- 
vatory (OVRO) millimeter interferometer and the nine 
6 m antennas of the Berkeley-IUinois-Maryland Associ- 
ation (BIMA) interferometer. The entire optical disk 
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of NGC 3521 (ll.O'xS.l') was mosaiced in 19 point- 
ings with Nyquist sampling of the 10 m antenna beam 
(FWHM of 1' for the 115 GHz CO J=l-0 line). The 
data were reduced and calibrated using the Multichan- 
nel Image Reconstruction, Image Analysis, and Display 
(MIRIAD) software package (|Sauh et al . 1995). 

We also obtained total power and short-spacing data 
with the 25-Beam Array Receiver System (BEARS) 
on the Nobeyama Radio Observatory 45 m telescope 
(NR045, FW HM = 15") Using the On-The-Fly observ- 
ing mode (Sa wada et al.|[200 8). the data were oversam- 
pled on a 5" lattice and then re-gridded with a spheroidal 
smoothing function, resulting in a final resolution of 
18.4". Data reduction was performed using the NOS- 
TAR data reduction package developed at the Nobeyama 
observatory. We constructed visibilities by deconvolving 
the NR045 maps with the beam function (i.e., a convo- 
lution of the 15" Gaussian and spheroidal function), and 
Fourier transforming them to the wy-space. After that, 
we combined the CARMA and NR045 data in Fourier 
space, inverted the uv data using theoretical noise and 
natural weighting, and CLEANed the maps. 

Throughout this paper, we adopt a CO-to-H2 con- 
yersion factor of Xq o=2.8xl0^° cm~^ (K km s^^)~^ 
(jBloemen et all [1989) for NGC 3521 and M51a. The 
validity of this Galactic factor is justified by the super- 
solar abundance averaged over the galactic disks (12 + 
log(0/ H) = 9.01 ±0.02 for NG C 3521, and 9.18±0.01 for 
M51a. IMoustakas et all 120101 ). Additionally, the small 
metallicity gradients across the disk justifi es the applica- 
tion of a single value (jBresolin et al.ll2Q0^ . 

For NGC 3521 at the distance of 8.03 Mpc, the angular 
resolution of the combined data (5.4"x4.7") corresponds 
to a physical scale of '^220 pc, and the r.m.s. sensitiv- 
ity (32 mJy beam~^) in the 5.08 km s~^ wide channels 
translates to a 3-ct level of 2.4 M0 pc^^ at this resolution 
(assuming i=72.7°). 

As for M51a, the detailed information of the CARMA 
-|- NR045 obs ervations and dat a reduct i on is described 
in iKoda et all (2009a) and Ko da et al.l ()20lH ). At the 
final resolution of the combined map (3.68" x 2.87", cor- 
responding to ~140 pc at a distance of 8.00 Mpc), the 
r.m.s. sensitivity (40 mJy beam~^) in the 5.08 km s^^ 
wide channels translates to a 3-a level of 18 M© pc~^ 
(assuming i=42°). For comparison, the typical size of 
a GMC in t he Galaxy is 50 to several hundred parsecs 
(jBlitz II1993D . We show a comparison between our com- 
bined CARMA + NR045 and the BIMA SONG maps 
at a matched resolution, 5.8" (corresponding to 230 pc 
for M51a) in Figure HJ This highlights both the higher 
depth of our image and the fact that we recover substan- 
tial faint-level emission, which is absent from the BIMA 
SONG image. 

2.2. Ha Emission-Line Images 

NGC 3521 was imaged using the 2.1 m telescope at 
Kitt Peak National Observatory as part of the SINGS 
ancillary data program (Ke nnic utt et al. 2003). The Ha 
line is contaminated by the doublet [N 11] A6548, 6584. 
A ratio of A6548, 6584/Ha=0.55±0.03 is adopted here 
to correct the data, following the integrated spectropho- 
tometric survey of nearby star-forming galaxies accom- 
plished bv .Moustakas et al.. (,201(^) . Our measurements 



on the final emission-line-only Ha image has a PSF with 
a 1.4" FWHM. The sensitivity limit of the emission-line 
image is 2.5 x 10" 



'^^ erg s ^ cm ^ arcsec ^ 



2.3. Spitzer MIPS Images 

Spitzer MIPS 24 fim maps are available for NGC 
3521, through the high level data products of the SINGS 
Legacy project (Spitz er Infrared Nearby Galaxy Survey, 
iKennicutt et al.ll2003D . The SINGS observation strategy, 
data reduction procedures, and map s ensiti vity limits 
are de scribed in IKennicutt et afl (|2003[ ) and IDale et all 
(|2005l ). For MIPS 24 ^m maps, the diffraction- limited 
angular resolution is 5.7" and the 1-cr sensitivity limit 
1.1x10^^ Jy arcsec^^. 

2.4. VLA H I Data 

The atomic hydrogen in NGC 3521 wa s mapped by 
the H I Nearby Galaxies Survey (THINGS, iWalter et all 
120081 ) using the NRAO Very Large Array (VLA). The HI 
maps for our study are obtained through the THINGS 
public data release websit^l. The H i map, as a re- 
sult of robust weighting, has a resolution of 8.19"x6.41". 
In each of the 5.2 km s~^ channels, the l-tr sensitivity 
is 0.47 mJy beam"'^, corresponding to a column den- 
sity of 5.8xl0^^cm~^ or 0.5 MqPC"^ in the integrated 
map. The calibration error for the H i data are esti- 
m ated to be ^10%. F urther technical details are given 
in IWalter et al.l ([20081) . 

3. METHODOLOGY OF SPATIALLY-RESOLVED S-K LAW 
STUDIES 

As mentioned in the Introduction, K07 and B08 find 
different results for the power law indices of the molecu- 
lar S-K Law in M51a. B08 qualitatively point out that 
the difference arises because K07 subtract a local Ha 
and 24 /im background for each aperture, but different 
photometric approaches (K07: aperture photometry vs. 
B08: pixel-by-pixel) may also play a role. Meanwhile, 
such a test involves a highly tricky issue of local back- 
ground removal. In this section, we present a new ap- 
proach of local background determination. By using it, 
we quantitatively show that removing or preserving the 
local background is sufficient to account for the discrep- 
ancy between K07 and B08. From here onward, we will 
indicate the fits to the Ssfr — relation as 



SsFR 



Mq yr~i kpc" 



= A 



M 



PC" 



7H2 



(1) 



with indicating the power law exponent of the molec- 
ular gas surface density. The dispersion about the fitted 
relation will be indicated with . The analogous rela- 
tion for the total hydrogen gas will be indicated as 



^SFR 



Mq yr~i kpc" 



oc 



M 



PC" 



(2) 



where E 



H 



Ehi + Eh, 



We also compare 



SsFR(FUV+24Mm) with EsFR(Ha-K24Ami) through the rela- 



tion 



SsFR(FUV+24^m) 0^ (EsFR(HQ+24Mm) 



\7SFR 



(3) 



^ http:/ /www. mpia-hd.mpg.de/THINGS/Overview. html 
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3.1. Local Background Subtraction 

In the hydrogen hne, FUV and mid-IR images used to 
trace SFR, besides the compact clumpy structures which 
trace of H ii regions and young massive stars, there ex- 
ists an underlying diffuse component of stehar/dust emis- 
sion, which is conventionahy referred to as "local back- 
ground" with slight ambiguity. This convention contin- 
ues in this paper, but we stress that the so-called "local 
background" here, apart from the local variation of sky 
or instrumental background, refers to the physical detec- 
tion of the emission from stars or interstellar dust. 

Pixel-by-pixel analysis introduces significant contami- 
nation from this local background especially in lower res- 
olution images, in which case the low luminosity pixels 
can be more strongly affected than those at the bright 
end, effectively overestimating the SFRs at the faint end 
and thus flattening the slope of the scaling laws. 

Conventional approaches to determine the local back- 
ground levels in galaxies include performing statistics in 
an annular region, or fitting the galactic image with an 
assumed two-dimensional analytical function (e.g., poly- 
nomial, exponential, Sersic profile or a combination of 
several galactic components). These methods work rea- 
sonably well in general but often fail in nearby galax- 
ies because of the large dynamical range and crowding 
of clustered regions. |Calzetti ct al. (2005) and K07 de- 
fined 12 rectangular regions in the M51a images that 
they used as their "background" areas. Since it is not 
sensitive to local variations on smaller scales, this ad-hoc 
strategy produces reasonable results (see next subsec- 
tion). However, a more automated approach is expected 
to enha nce the repr oduceability and robustness of the 
results. I Blanc etall (pOOS) identify and remove the dif- 
fuse ionized gas (D IG) from the Ha ima ge of M51a sta- 
tistically, based on iMadsen et al.l ()2006[ ) which demon- 
strate that the DIG in the Galaxy has a differenct [S 
II] /Ha ratio (0.34) from that of Galactic H ii regions 
(0.11). The iBlanc et al.l (j2009) procedure has the ad- 
vantage of being astrophysically motivated, but the DIG 
fraction recovered by these authors (11%) is smaller than 
that recovered by photometric m easurements either when 
the whole galaxy (40-50%, e.g.. iGreenawalt et allll998l: 
iThilker eF al. 2000) or the central part of M51a is consid- 
ered (we find the f raction to be 32 % in the same central 
part covered by Blanc et al.| l2009l . using Hllphot). 

In this work, we distinguish the compact component 
(presumably H ii regions) from diffuse stellar/dust emis- 
sion through morphological examination. This strategy 
is inevitably arbitrary to some extent, but its basic moti- 
vation has a strong physical underpinning. The proper- 
ties of the diffuse Ha, 24 /im and UV emission have been 
extensively discussed by the existing literature, which is 
briefly reviewed in the discussion section. When such a 
strategy is being applied, the crowding of sources is the 
major obstacle to accurate background determination, 
but it can be circumvented by removing the compact 
clumps instead, which is done by adapting the algorithm 
of the IDL software Hllphot developed by IThilker et al.l 
((2000h for our purposes in this work. The algorithm can 
be summarized as follows. On an Ha image of a galaxy, 
it identifies H ii regions by smoothing the image with 
kernels of different sizes, and utilizes object recognition 
techniques to identify significant peaks ("seeds"). After 



that, an iterative procedure is started allowing for area 
growth from these "seeds" . The boundaries are built up 
at successively fainter isophotal levels that keep grow- 
ing until either boundaries of other regions are encoun- 
tered, or a pre-established lower limit to the gradient in 
surface brightness is reached implying its arrival at the 
background level. Two-dimension interpolations are then 
carried out by fitting the background pixels around each 
H II region to estimate the local background behind the 
source. 

Originally designed for automated photometric char- 
acterization and statistical analysis of H ii regions and 
their luminosity functions in Ha images, the algorithm 
of Hllphot can be easily adapted to 24 ^m and FUV 
images. We tested our implementation of HHphot on 
the HST Ha image of M51a, and produced an H ii re- 
gion l uminosity func ti on which is very similar to that of 
.Scoville et all (|200lD ()Liu et al.l 120 ll( ). It should be no- 
ticed that our purpose of extracting star- forming regions 
against the local background is somewhat different from 
conventional Ha data analysis which aims at removing 
any contamination from diffuse components in H ii re- 
gions. The lower limit gradient in surface brightness for 
the iterative growing procedure from the "seeds" is thus 
chosen to be sufficiently low (lower than what is typically 
chosen for H ii region analysis) that we are confident to 
have reached an appropriate background level, but still 
sufficiently above zero to prevent growing deep into am- 
bient areas. 

Because of the much lower contribution from the dif- 
fuse emission in Ha than in mid-IR and FUV, the re- 
moval of the Ha local background is done conservatively, 
aiming at a background level to be reached only in the 
ambient area between the outer spiral arms. For the 
purpose of testing this approach, we also use a constant 
global background determined using several external H 
II regions in deriving S-K laws, and the resultant slopes 
differ minimally (< 0.1, typically ^^0.05). 

The established source-free image is not smooth, but 
contains sharp edges at the boundaries of sources. We 
therefore smooth the image by calculating the median of 
the pixels enclosed in a moving box and replacing the 
central pixel with the median. Because of the robust- 
ness of the median in statistics, this smoothing technique 
leads to a relatively stable result and suffers much less 
from the steep local variations which will affect an im- 
age smoothed by kernel convolution. The size of the 
box is chosen as a compromise between minimization 
of the crater-like artifacts left by the removed H ii re- 
gions and minimization of detail loss. The final result- 
ing "diffuse emission" map is then ready for subtraction 
from the original image. The background-free images, 
along with the background images before and after the 
median smoothing, are shown in Figure [3] for the pur- 
pose of comparison. As can be seen from the figure, the 
background-free image mostly contains clustered sources, 
and the disk-like diffuse emission has been removed. 

3.2. Linear Relation without Local Background 
Subtraction 

BOS performed a pixel-by-pixel analysis on 7 spiral 
galaxies detected in CO (J=2-l) at a resolution of 750 
pc including M51a. The employed SFR tracer in their 
study is a combination of far-UV and 24 /xm luminosities 
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without local background subtraction and whose coef- 
ficients are adjusted to be consistent with the Ha-|-24 
/xm SFR tracer at 750 pc resolution, th e applicability 
of wh ich has been justified in K07 and iCalzetti et alJ 
([2007'). No correlation between SFR and H i surface 
density is seen, a similar result to K07. In contrary to 
the super-linear relationship found by K07 (power law 
index = 1.37 ± 0.03), BOS find the surface den- 
sity of SFRs to be proportional to that of molecular gas 
(power index 7H2 = 1.0 ± 0.2). We notice that BOS car- 
ried out their analysis using CO (J=2-l) images from 
The HERA CO L ine Extragalactic Survey (HERACLES; 
iLerov et al.ll2009l) except for M51a, for which they share 
the same data set as K07, including the CO (J=l-0) 
image from BIMA Surve y of Nearby Galaxies (BIMA 
SONG: lHelfer et al.l[200l . This enables us to make a di- 
rect comparison of the two approaches of measurements 
adopted in these two studies to identify the reason for 
the discrepancy of the obtained power law indices. 

The KPNO Ha, SINGS 24 fim, GALEX FUV, 
THINGS 21-cm H i and BIMA SONG CO (J=l-0) 
images of M51a have resolutions of 1.9", 5.7", 4.6", 
5.S2"x5.56" and 5.S"x5.1", respectively. The 3-ct sen- 
sitivity is 1.3 Mq pc~^ for the H I map (K0 7), and 13 
Mq pc-2 for the CO map (jHelfer et al.l[200l . The 1-cr 
detection limits of the 24 /zm (l.lxlO~^ Jy arcsec"^). 
Ha (l.SxlO"^"^ erg s~^ cm~^ arcsec"^) and FUV images 
(3.6x l0~^^ erg s~^ cm~^A~^ arcsec"^) (jCalzetti et al.l 
I2OO5I ) lead to a corresponding SFRHQ+24^m of 1.0x10 
Mq yr"i kpc~^ and a SFRFuv+24A.m of 1.7 x 10"^ 
Mq yr~-^ kpc~^ at 750 pc (19 .3") reso l ution Here we 
have assumed an IMF given bv lKroupal ()200lD to derive 
the SFR surface densit y (in units of Mq yr~^ kpc~^) 
following Calzetti et al.' (j2007[ ) and ILerov et al.l (HOOS): 

SFRHa+24Mm = 5.3 X 10-^^ (^^^^ ^ o_o31 L24pm), (4) 
SFRFUV+24pm = 3.4 X 10-4^ (ipUV + 6.0 i24Mm), (5) 

where i24/^m and Lpuv, both defined as XL\ and in units 
of erg s^^, are referenced to the central wavelength of 
the Spitzer 24 ^m (A=23.6S^m) and the GALEX FUV 
(A=152SA) filters, respectively. 

We first convolve these M51a images with circular or 
elliptical 2-d gaussian kernels so that all the images are 
downgraded to a common physical resolution of 750 pc, 
the working scale of BOS, for the purpose of comparison. 
Aligned to the astrometric frame of the 24 fim image 
which has the coarsest grid (1.5"), they are then regrided 
to a pixel scale of 750 pc to avoid sub-PSF sampling. 
The data are then converted to physical quantities, cor- 
related pixel-by-pixel, and analyzed by linear fitting. In 
this study, the linear regressions are performed with the 
ordinary least square bise ctor method (same as K07 and 
BOS), as recommended bv llsobe et al.l ()1990f ). 

The results for M51a are summarized in Figure |4] (the 
left of the panel pairs). At this resolution, we reproduce 
the resuhs of BOS (7h2 308=O.S4) when FUV + 24 ^m 
SFRs are correlated to molecular gas surface densities 
(7H2=0.S3±0.04, Figure IH-b). we find a similar sub- 
linear correlation (7h2=0.86±0.05) when using Ha-|-24 
fim as a SFR indicator for the molecular-only S-K law 
(Figure ID-a), because of the tight, close-to-linear corre- 



lation between SFRFuv-i-24;im and SFRHQ+24;im in the 
SFR > 10~ ^ Mq yr~^ ^^P'^"^. regime (a result of the cal- 
ibration by |Lero^^F^L| [2001, Figure|4]-c). Nevertheless, 
the two SFR tracers correlate poorly with each other 
when SFR is lower than 10"'^ Mq yr~^ kpc~^, hinting 
at a discrepancy of deriving Ha and UV at low surface 
brightnesses. All these linear fits have taken into account 
only those data points above their corresponding 3-ct sen- 
sitivity limit which are shown in each panel as horizontal 
or vertical dotted lines. 

Little or no correlation is found between SFRFuv+24;im 
or SFRHQ+24;^m with the atomic gas surface density, and 
we observe a saturation at Ehi '~ 10 Mq pc~^. Both 
effects had been already observed by both K07 and BOS. 
An evident "kink" is seen at the transition from atomic- 
to molecular-dominated regimes when both components 
are combined to make a SFR - Sh comparison (Figure H]- 
e andlD-f). Compared to Figure |4]-a or |4]-b, the slope 
hardly changes in Figure [H-e (7h=0.S6±0.03, dashed 
line) where we restrict the fitting to pixels in which the 
total hydrogen surface density and SFR are above 3ct. 
However, when extended to the outer disk by assum- 
ing little or no existing molecules, the linear fit gives 
7H=1.26±0.02, steeper than 7h=1.11 found by BOS, a 
consequence of the fact that we do not mask the regions 
outside D25- The results for M51a presented by BOS 
are therefore thouroughly reproduced by our indepen- 
dent pixel-by-pixel analysis. 

We perform a similar comparison to BOS's result for 
NGC 3521, ahhough we use our CARMA + NR045 
CO ( J=l-0) map, while those authors used their HER- 
ACLES CO (J=2-l) data. Since we adopt a smaller 
distance for this galaxy than BOS (8.03 Mpc vs. 10.7 
Mpc for those authors), their assumed 750 pc pixels cor- 
respond in reality to 1 kpc pixels. We thus re-grid NGC 
3521 to 1 kpc pixels and obtain a molecular-only power 
law index 7H2 = 0.90 ± 0.13 when SFR is derived from 
Ha and 24 ^m data, or 7H2 = 0.92 ± 0.12 when FUV 
is used instead (BOS find 7H2308 = 0.95). The SFRs 
derived in these two different ways coincides with each 
other, SFRFuv+24pm « SFRjj^™m- The scaling rela- 
tion for total hydrogen gas content gives 7h — 2.06±0.31 
(without masking any part of the H i map), similar to 
BOS's value of 2.12. We infer that the two different CO 
maps (J=l-0 and J=2-l) do not lead to discrepancies 
in the resulting S-K law, at least at the spatial scales 
investigated for this galaxy. 

3.3. Super-Linear Relation with Local Background 
Subtraction 

By employing both the extinction-corrected Paa line 
image and the weighted sum of Ha and 24 fim luminosi- 
ties as an unbiased SFR tracer, K07 performed photome- 
try using apertures with 13" diameter {— 520 pc) on 257 
positions in M51a where star formation is detected. The 
authors identified 12 rectangular areas and fit the local 
background for subtraction. This strategy was applied to 
both the hydrogen line image and the 24 /im map. As a 
result, a best-fit slope of 7h = 1.56±0.04 was found when 
the specific SFR is correlated to the total gas surface den- 
sity on 520 pc scales; the slope decreases to 1.37±0.03 on 
1850 pc scales. Additionally, virtually no correlation is 
seen between the local SFR surface density and the H i 
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surface density which saturates at log [Ehi/M© pc~^]~l, 
while a strong correlation with a shallower best-fit slope 
of 7H2 = 1-37 ± 0.03 exists for the molecular S-K law on 
520 pc. 

The actual consistency of the results from K07 and BOS 
can be shown using the background-subtracted 24 /im, 
FUV and Ha images of M51a, combined with the BIMA 
SONG CO map. We switch our pixel-by-pixel analy- 
sis to 507 pc (adopting a distance of D=8.0 Mpc), the 
scale examined in K07 (520 pc, assuming D—8.2 Mpc), 
to retrieve their results, where OLS bisector fittings are 
performed again, with the local background subtracted 
in the SFR map. Our 500 pc molecular-only S-K law in 
M51a has a slope of 7H2 = 1-48 ± 0.06, slightly steeper 
than but only 1.3(7 away from what is obtained by K07 
through Ha + 24 ^m SFR, 7h2,ko7 = 1-37 ± 0.03. The 
best-fit power relation between the two SFR tracers is 
given by SFRfuv+24 « SFR^f+±4^°^. This calibra- 
tion now has been improved relative to that from the 
unsubtracted images (section 3.2) by a remarkably weak- 
ened scatter at the faint end, and the tight correlation 
now extends down to SFR ~ 10""* M© yr~^ kpc^^, im- 
plying the applicability of FUV -|- 24/xm as a SFR tracer 
as long as the local background is properly removed. Not 
surprisingly, changes insignificantly to 1.47±0.05 if 
unobscured SFRs are measured by substituting Ha with 
FUV luminosity. A even stronger sign of no correlation 
between Ehi and Ssfr is seen. For the total gas, we 
find a best-fit 7h=1.52±0.05 when a 3-cr sensitivity limit 
cutoff is adopted for the total hydrogen mass surface den- 
sity, virtually identical to the K07 result. However, from 
a practical point of view, only those data points with 
both H I and CO detections above their respective 3a 
should be taken into account for the linear fit, just as 
K07 did in their study of total gas surface density scal- 
ing relation. Such a criteria results in a steeper slope 
7h=1.66±0.06, still consistent with K07. Hence, all the 
results that K07 obtained by performing aperture pho- 
tometry are reproduced with our pixel-by-pixel analysis 
on the background subtracted images. 

In conclusion, using the same set of data as K07 and 
BOS, our analysis accurately reproduces the results of 
both investigations. We hereby conclude that at least for 
the case of M51a, the crux of their discrepancy lies in the 
preservation/elimination of the local background in 24 
/im, FUV maps and Ha images, ruling out the different 
employed measuring strategies (aperture photometry vs. 
pixel-by-pixel analysis) and any other factors as possible 
important causes. 

4. STAR FORMATION IN M51A AND NGC 3521 ON 
SUB-KPC SCALES 

The BIMA SONG CO map of M51a has been used so 
far for the purpose of comparison only, and is now re- 
placed with the more recent CARMA + NR045 data to 
actually probe the spatially-resolved S-K law. For NGC 
3521, the spatial resolution of our data set is limited by 
that of the H i map (S.19"x6.41"), corresponding to a 
physical scale of 320 pc at the adopted distance, S.03 
Mpc. However, the resolution limit of molecular-only S- 
K law can be pushed to smaller scale values, to 220 pc or 
5.7", the PSF FWHM of the 24 /xm image. As for M51a, 
the similar resolutions of its H i map (5.S2"x5.56") and 
24 fj,m image set the accessible limit to be ~ 230 pc. 



The H I and 24 /im maps, which have the coarsest pixels 
(1.5"), define the grid onto which images are registered. 
We then smooth all the images to make the PSF circular 
with the same resolution and rescale the pixels to match 
that resolution before converting flux and surface bright- 
ness units to physical quantities. An inclination angle 
of 72.7° (42°) results in a correction factor of 0.29 (0.74) 
for projected surface densities in NGC 3521 (M51a). The 
CO data, having a lower sensitivity limit than H i by a 
factor of ~2.3, set our limiting sensitivity to the cold gas 
surface density measurement. Justified by the results of 
the previous sections, we now apply a pixel-by-pixel anal- 
ysis to both galaxies, M51a and NGC 3521, utilizing SFR 
images where the diffuse emission/local background has 
been removed. The S-K law of M51a and NGC 3521 are 
shown in Figure [5] and Figure [5] at their highest available 
resolution, respectively. 

We now investigate the scaling relationships at differ- 
ent spatial resolutions. The dependence of best-fit pa- 
rameters on spatial resolution is shown in Figure[7l where 
a series of power-law indices for M51a and NGC 3521 are 
shown in parallel. The range of spatial scales is different 
for these two objects due to their difference in inclination 
which yields a factor 2.5 smaller apparent minor axis in 
NGC 3521 (^3 kpc across) than in M51a. Investigations 
on scales larger than a quarter of this length will be of lit- 
tle significance because of the paucity of data points and 
of including a too large fraction of the galaxy. Moreover, 
the high inclination of NGC 3521 makes the removal of 
the background non-star-forming light more complicated, 
which can affect our results especially on relatively larger 
scales. Considering the inclination of NGC 3521, 700 pc 
actually corresponds to a physical scale of ~2 kpc, which 
we take as a convenient limit for our study. The best-fit 
parameters of the molecular S-K law in M51a and NGC 
3521 at different resolutions are reported in Table [Hand 
Table [5J where SFR is derived from the weighted sum of 
Ha and 24 /im luminosities. 

Figure [7] demonstrates that 7sfr has little dependence 
on spatial resolution in both galaxies, varying between 
0.97 and 1.01 in M51a and between 0.99 and 1.04 in 
NGC 3521. The weak ascending trend of 7sfr is caused 
by the larger scatter at the faint end as discussed in the 
methodology section, and the uncertainty in background 
removal may also play a role. The FUV-l-24 /im lumi- 
nosity is therefore also a good SFR tracer as long as the 
diffuse emission component is properly removed. 

Figure[7]also reveals that the power-law indices of both 
the molecular and the total hydrogen S-K law decrease 
monotonically as the spatial scale increases. Specifically, 
as the resolution decreases from 230 pc to 1 kpc, 7H2 
decreases from 1.9 to 1.2, implying higher star formation 
efficiency (defined as the SFR-to-gas surface density ra- 
tio) in the denser regime of molecular gas on all sub-kpc 
but super-CMC scales. For the total gas, we find 7h de- 
creases from 2.1 to 1.6 on the whole range of scales. At 
'--^500 pc resolution, wc find steeper slopes (7H2 ~ 1-5, 
7h ~ l.S) than K07 (7h2,ko7 ~ 1-4, 7h,ko7 ^ 1-6), 
mainly because the new CARMA -I- NR045 data is more 
sensitive than the BIMA SONG map by a factor of ~2, 
and has collected more diffuse emission from CO, push- 
ing the low density end to higher values. A even steeper 
slope would be found if we followed BOS, assumed no 
molecular gas outside the star-forming disk and fit the 
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scaling relation on the whole isophotal area defined by 

025. 

In contrast to M51a, the high inclination of NGC 3521 
results in blending of star formation sites which compli- 
cates the removal of diffuse stellar/dust emission, and 
can affect our result on relatively larger scales. Figure [7] 
shows that the slopes of the S-K law of both the total 
and molecular gas and of the SFR indicator comparison 
for NGC 3521 are less certain and stable than M51a, 
indicating higher uncertainties in the local background 
subtraction. In addition, the FoV of our CARMA obser- 
vation covers only about half of the D25 scale, missing 
the information from the outskirt of the disk where the 
transition from dominant molecular phase to dominant 
atomic phase in the cold gas occurs. This shortcoming 
tends to reduce the amount of data points in the regime 
slightly denser than Sh '~ 10 Mq yr"-'^ kpc~^ and could 
have steepened the resulting slope. Bearing in mind these 
potential issues, we still observe a descending trend sim- 
ilar to that of M51a in both 7H2 and an^ . However, both 
quantities are smaller than in M51a at a given resolution: 
7H2 is flatter by ~ 0.3, and cthj smaller by ~ 0.15. The 
power-law index of 7h in NGC 3521 fluctuates weakly 
within the errors and persists between 2.1 and 2.2, im- 
plying very weak or no dependence on the resolution. 
However, as we have mentioned above, this may due to 
the absence of the CO data in the outer star-forming disk 
which hinders a complete sampling over the range of gas 
surface density. 

The scatter or dispersion of the S-K law, quantitatively 
defined as the r.m.s. of the fitting residuals measured 
perpendicular to the best-fit line, is an important param- 
eter indicating the validity of the S-K law. The study of 
K07 suggests that the 500-pc-scale total-gas S-K law has 
a scatter almost twice as large of the 2-kpc-scale mea- 
surements in M51a, which becomes the first hint of the 
existence of some scale below which the total-gas S-K law 
breaks down. Figure [7]-b reveals the monotonic decreas- 
ing trend of the scatter ((TH2) about the best-fit molec- 
ular S-K law as the resolution becomes lower on sub- 
kpc scales in both our program galaxies. iMomose et al.l 
mm) analysed the CARMA + NR045 CO data of NGC 
4303 and find that the K-S law already breaks down at 
~250 pc resolution, and does not recover unless the data 
are smoothed to a scale of 500 pc or larger, especially in 
the disk region of the galaxy. Converselly, iVerlev et al.l 
((2010:) show that a loosely correlated tota l S-K law exists 
down to 180 pc in M33 (0.84 Mpc), and lOnodera et al.l 
(|2010f) conclude that the S-K law of H2 in M33 becomes 
invalid only at the scale of GMCs ('-^80 pc). Limited by 
the angular resolution, we are unable to reach scales be- 
low ~200 pc, but significant scatter (crH2=0.62 dex for 
M51a and 0.52 dex for NGC 3521) about the best-fit 
molecular S-K law is observed in both galaxies at their 
highest available resolution. When the linear size in- 
creases by a factor of 3 (^700 pc), the scatter decreases 
to 0.48 (M51a) and 0.26 dex (NGC 3521). Assuming the 
trend in Figure [7] [b] extends to smaller scales, a rough 
extrapolation gives a scatter of close to 0.7 dex at 80 
pc scale, enough to mask any trend. The smoothly in- 
creasing scatter suggests the absence of a characteristic 
breakdown scale of the S-K law, but the scaling rela- 
tionship will be overcome by the large scatter gradually. 
Higher resolution datasets are required to definitively pin 



down this issue. 

As an exploratory experiment, we also attempt to find 
a common formulation that relates the observational re- 
sults in the two galaxies, so that some hints for the un- 
derlying physical universality may be obtained. Since 
the molecular S-K law is what should be focused on on 
sub-galactic scales, 7H2 is of interest in the first place. 
Furthermore, because the uncertainty in the photometric 
measurements is typically at the level of ~10-20%, the 
dispersion about the best-fit S-K law is mainly intrin- 
sic, as has been pointed out by previous studies (K07; 
iBlanc et al.ll2009t ). and 

(JH2 is a quantity of physical im- 
portance. 

The uncertainties of the slopes 7H2 in Tables [1] and [2] 
are, in fact, fitting errors, and are consistent i n value with 
those already reported in the literature (K07: lBlanc et all 
l2009f l. These uncertainties are not fully related to the 
dispersions tTH2i since, for instance, a larger dynamical 
range on the data could decrease the slope uncertainty, 
even at constant cth2 ■ Hence the relevance of the intrinsic 
data scatter an^ , which we report as separate columns 
in Tables [Hand [2 

As mentioned above, 7H2 and au.^ are systematically 
smaller in NGC 3521 than in M51a. However, the reso- 
lution under consideration is the projected scale, and the 
comparison between the two galaxies should be made on 
the basis of physical or de-projected scales. We hence re- 
plot the dependence of 7H2 and o'h2 on Sdp (= <5/cos i) 
in Figure m As revealed by the two panels, on sub-kpc 
but super-CMC scales, both quantities are surprisingly 
consistent between the two galaxies. This remarkable re- 
sult is a sign that we are approaching the scale where 
the physical origin of the S-K law is appearing. We are 
motivated to propose the following relations which hold 
in both galaxies: 

7H2 = A log [Jdp/kpc] + B, (6) 

and 

tTH2 = C [(5dp/kpc] + D, (7) 

where cth2 is in dex, and A — —1.09 ± 0.05, B — 
1.36 ± 0.02, C = -0.19 ± 0.01 and D = 0.67 ± 0.01. 
These best-fit parameters are determined using the data 
of M51a only, which are superior than those of NGC 3521 
in many respects, as explained above. The values change 
minimally if the latter is included in the fits because of 
the impressive overlap among the data (Figure [8]). This 
overlap seemingly indicates a universal sub-kpc molecu- 
lar S-K law that persists amongst spiral galaxies whose 
actual form depends only on the physical scale being con- 
sidered. An extrapolation of the above expression for 7H2 
shows that a slope of unity requires a scale as large as 
2 kpc, roughly the typical spacing of adjacent spiral 
arms in disk galaxies similar to M51a and NGC 3521. 
This has been verified in NGC 3521 (Figure [7]-a), where 
7H2 is consistent with unity within the error bars at the 
apparent scale S ^ 700 pc (Tabled]), which corresponds 
to a de-projected scale Sdp ^ 2 kpc. 

Similarly, we find that the - S correlation in M51a 
can be fitted to 

7H = -0.88 log [S/kpc] + 1.60. 

For a projected scale (5=1.80 kpc (corresponding to 1.85 
kpc if D=8.2 Mpc is adopted as in K07), the extrapola- 
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tion of the above function predicts 7h=1-36, well consis- 
tent with the slope 1.37 ± 0.03 found at this resolution 
by K07 for the total hydrogen S-K law using the 14-m 
Five Col lege Radio Astronom ical Observatory CO ( J=l- 
0) map ([Lord &: Younel ll990l. Despite this consistency, 
we are aware of the uncertainties in deriving large scale 
information from the extrapolation of our sub-kpc study. 

The discussion above leads to a conclusion that caution 
must be taken in interpreting observational results, es- 
pecially those of galaxies with high inclinations, which 
makes the intended sub-kpc studies actually "super- 
kpc". Linear or sub-linear molecular S-K laws may be 
found if spiral structures are not readily resolved. The 
HERACLES CO (J =2-l) data that i s used by BOS has 
a resolution of 11" (jLerov et al.l 12009? ) . This spatial res- 
olution (430 pc at a distance of 8.03 Mpc where NGC 
3521 is located, or 570 pc at 10.7 Mpc adopted by BOS) 
translates to Sdp ~ 1.5 kpc, already close to the ~ 2 
kpc unity-slope scale, and is virtually incapable of re- 
solving the spiral structure of NGC 3521, which is likely 
to be another important reason for the results in BOS in 
addition to preserving the local background in the SFR 
measurements. 

At our highest available physical resolution ((5dp=320 
pc) in M51a, the power indices of both the molecu- 
lar and total hydrogen S-K laws are about 2 (7H2 — 
1.91 ± 0.03, 7H = 2.14 ± 0.04). A slope as steep has 
never been reached in spatially-resolved studies of ex- 
tragalactic star-forming disks, but has been observed 
inside the Milky Way. For instance, iMisiriotis et al.l 
(|2006[) derive both the gas/dust distribution and SFR 
from COBE/DIRBE da ta and find 7h = 2. IS ± 0.20; 
iGutermuth et al.l ()2011f ) directly observed about 7000 
young stellar objects (YSOs) and obtain a slope ^2, 
where the gas density is traced by dust extinction. Our 
results in M51a are therefore very similar to these Galac- 
tic investigations, which seems to hint at an "intrinsic" 
S-K law with slope ^2 in spiral galaxies, and the flat- 
ter power indices that have been found in other galax- 
ies are actually a resolution effect. As a minimum, we 
have added evidence that a super-linear relation is found 
in extragalactic investigations supporting the recent re- 
sults that the star formation efficiency positively cor- 
relates with the column densit y for Galactic molecular 
clouds ()Gutermuth et al.ll2011|) . As a caveat, our study 
has been performed only on two galaxies, and more data 
with better sensitivity, resolution and broader FoV for 
a larger sample of galaxies are required to confirm the 
results presented in this section. 

5. DISCUSSION 

5.1. The Nature of the Diffuse Emission at all 
Wavelengths 

Among the SFR images, the extended background 
contamination is th e most significant at 24 /im. 
iPopescu et~aLl ()2005D find that the IR-to-UV ratio is 
higher in the interarm regions of MlOl than the arm 
regions, which they interpret as due to dust heated by 
optical photons in the interstellar radiation field. In 
order to account for the observed SEDs of dust emis- 
sion, a dust model must include a substantial population 
of ultrasmall grains or large molecules with the vibra- 
tional properties of PAH material and with sizes such 



that single-phot on heating can exci te the observed vibra- 
tional emission (|Draine "fcLil[2007l) . Therefore, the micl- 
IR emission is produced not only by dust grains heated 
by young stellar populations, but also, and in some areas 
predominantly, by dust heate d by o lder stars through 
single-photon processes fe.g.. IDraine fc Li 2007). The 
latter component is unassociated with current star for- 
mation, and is likely to dominate the mid-IR emission in 
the faint, diffuse regions in galaxies. As a quantification 
of the amount of diffuse light at each position in the two 
galaxies, we show in Figure |9] the ratio of the compact 
(background-subtracted) emission to the diffuse compo- 
nent (the subtracted background) in Ha, 24 /xm and FUV 
as a function of H2 mass surface density at 400 pc reso- 
lution (projected scale). It can be seen that, among the 
three wavelengths, the 24 /im emission shows the most 
marked trend for increasing compact-to-diffuse ratio with 
increasing (especially in M51a). The variation of the 
compact-to-diffuse ratio is broader at 24 /im than Ha and 
FUV (especially in M51a), in line with the strongest in- 
fluence from its diffuse emission. In general, the fact that 
young stars form in dusty galactic environment leads of- 
ten to the prevalence of dust-obscured SFR (traced by 
IR) over unobscured SFR (traced by optical and UV), 
implying that the slope of the S-K law is mainly driven 
by the way the diffuse 24 /im emission is handled. 

In the FUV data, the diffuse compone nt is also con- 
sider a bly stronger than i n the Ha images (jMeurer et al.l 
119951 iMaoz et all 119961 ). likely owing to the longer 
timescale over which UV continuum photons can be pro- 
duced. As an example, it takes only 9 Myr for the Ha 
emission of an instantaneous-burst population to dimin- 
ish by two orders of magnitude, but takes about 100 
Myr at 1500 A where the FUV band of GALEX resides 
(jLeitherer et al.l I1999D . As a result of galaxy dynam- 
ics, over 100 Myr, the UV-emitting stars will have mi- 
grated away from their birth-site, and will have diffused 
across the galactic disk over an area at least 10 times 
larger than that of the ionizing-photon-cmitting stars 
(Chandar et al. 2005; Pellcrin et al. 2007). This results 
in a vastly reduced contrast between the arm and inter- 
arm regions of galaxies in FUV (~15:1 in M51a after ex- 
tinction correction) relative to the Ha (>30:1 in M51a), 
and also in a less obvious association between the loca- 
tion of the FUV emission and present-day star formation. 
Investigations on the properties of diffuse UV emission 
have demonstrated that, the diffuse UV light is likely 
originating from evolved, disrupted stell ar clusters that 
has migrated away from the ir birth-sites (jTremonti et al.l 
I2OOII : iChandar et al.l [20051 ): any assumptions of in-situ 
star formation for the diffuse UV light implies a steep 
IMF which is either steeper than Kroupa at the high 
mass end (>10-20 Mp) or t runcated (down to 30-40 M©) 
(e.g., iTremonti et al.ll200lD . 

To estimate the impact of the diffuse UV light, we as- 
sume that star formation proceeds at a roughly constant 
rate over 1 Gyr, and that stars and mol ecular clouds 
remai n closely associa ted over ~30 Myr (|Bhtz fc Shu I 
il980tiElmegreen ||2007[ ). in agreement with some observa- 
tions for the oldest stars found in GMCs ()01iveira et al.l 
|2009() . St ellar population synthesis models (STAR- 
BURST99. iLeitherer et~all [r999) indicate that for con- 
stant star formation, S0% of the GALEX FUV light is 
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produced by stars younger than 30 Myr, and 20% by 
stars older than 30 Myr. If the older stars get dis- 
persed over the galactic disk, their mean dust atten- 
uation will be lower than that of the younger stars, 
which are chistcrcd in the spiral arms. Using the opacity 
values of .Holwcr da et al.l (pbOSi ). A/^arm ^^1.5 mag and 
^/.intorarm ~0.5 mag, and the conservative assumption 
of a mixed dust /star geom etry, with a Milky Way ex - 
tinction curve plus albedo ()Calzetti et al.l 11994 I2OOO0 . 
the contribution to the total observed FUV emission of a 
galaxy from stellar populations younger than 30 Myr is 
only about 60%. Other assumptions for the dust/star ge- 
ometry and the effective attenuation will lead to smaller 
fractions contributed by the young star populations to 
the total observed FUV. In light of the above modeling, 
unless stellar populations significantly older than 30 Myr 
remain closely associated to their birth cloud, at least a 
40% fraction of "diffuse" UV light should be present. Our 
strategy of diffuse emission subtraction presented here re- 
sults in an observationally determined diffuse FUV frac- 
tion of 44% in M51a and 53% in NGC 3521. By virtue of 
its lower inclination, M51a has a more reliable measure- 
ment, and our modeling and observation are therefore 
consistent. 

Although removing the diffuse FUV emission is clearly 
a necessity for measuring SFRs, we have performed a 
test by making a SFR map of M51a with local back- 
ground preserved for the FUV image but subtracted for 
24 /im. Compared to our fiducial method of SFR trac- 
ing (IlQ;-|-24/L(m with the background removed in both 
bands), this SFR map reduces the slope of the molec- 
ular S-K law from 1.42±0.05 to 1.14±0.04 at 750 
pc resolution, higher than the slope 0.83±0.04 obtained 
when both the FUV and 24 /im backgrounds are pre- 
served. Thus, keeping the FUV background recovers a 
slope roughly half-way between those obtained when ei- 
ther removing or preserving both backgrounds. How- 
ever, we still recover non-linear trends at high spatial 
resolution: 7h2 = 1-39±0.02 at 250 pc resolution, when 
preserving the FUV background and removing it from 
the 24 fim image, to be compared with 7h2=1.86±0.03, 
when the background is removed from both images. 

The stronger diffuse emission in UV relative to the 
Ha is the underlying reason for the deviation of the 
SFRFuv+24/jm vs. SFRHa+24/im Correlation from a sin- 
gle power law at the faint end (see Figure |4]-c) . As we 
have shown, this deviation can be calibrated by careful 
local background removal, whose necessity is thus justi- 
fied by the need to include in our SFR indicators only 
the young stellar populations (and the dust they heat) 
involved in the current star formation activity. By in- 
cluding the local background contamination in their SFR 
estimates, BOS systematically overestimates the current 
SFRs at the faint end, which affects the physical inter- 
pretations of their high quality observations. 

The Ha data are processed with our strategy as well, 
driven by the fact that a fraction of ionizing photons es- 
cape from H II regions and result in the diffuse emission 
from ambient ionized gas unrelat ed to the current star 
form ation at a sp ecific location (iFerguson et al.l I1996L 
Il99a iWang et al.[|1999l: iHeckman et al.l 119990 . Another 
reason is that pixel-by-pixel analysis is highly sensitive to 
the low level uncertainty in flat-fielding, continuum sub- 
traction and sky background removal and requires careful 



handling. In fact, this technical difficulty can be circum- 
vented by u s ing in tegral field spectroscopy, as is done by 
I Blanc et al.l ()2009D . The authors investigated the central 
4x4 kpc^ of M51a at 170 pc resolution, with the diffuse 
Ha emission corrected based on [S ii]/Ha line ratio, as an 
alternative to the local background removal in K07 and 
our work. Combining these data and the archival BIMA 
SONG CO map, they fit the data through a Monte Carlo 
approach with upper limits included and find a sub-linear 
slope — 0.82 ± 0.05 for the molecular S-K law, sim- 
ilar to 0.84 obtained by B08. However, as they have 
pointed out, when an ordinary least-square (OLS) bisec- 
tor linear fitting excluding upper limits (as throughout 
this paper, K07 and B08) is performed on their data, the 
slope steepens to 1.5, close to ~1.4 found by K07. At 
our highest available resolution (230 pc) when the same 
BIMA SONG map is used, we find a similar power index 
1.45 ± 0.03 if we on ly exclude all data points with signif- 
icance below 1(7 as iBlanc et al.l (2009) did, rather than 
our default 3(T. Hence, they actually reach a conclusion 
consistent with K07 and this study, when noisy (below 
la) CO data are removed from their analysis. 

5.2. Robustness Tests 

In this work, we have performed our analysis in a cus- 
tomary manner, i.e. b y working above a certain sensi- 
tivity threshold (3-(t). IBlanc et al.l ()2009D simulate the 
likely distribution of non-detections so that data below 
1-fT detection limits can be included in their power-law 
fitting. While this approach attempts to avoid data cen- 
soring, their strategy requires extremely accurate under- 
standing and reliable modeling of the noise characteris- 
tics, which is a challenging enterprise especially for CO 
maps. Unfortunately, the current major limitation to 
spatially-resolved S-K law studies is the relatively low 
sensitivity of the CO maps, and the real solution hinges 
on the higher quality of future CO data (through, e.g. 
ALMA). Our CARMA + NR045 data for M51a is more 
sensitive than the BIMA SONG map that IBlanc et al.l 
(2009) used by a factor of ~2, and has higher image 
fidelity. Figure [2] compares these two maps at 230 pc 
scale (corresponding to the BIMA SONG resolution), 
and demonstrates that the CARMA -I- NR045 observa- 
tions recover slightly more flux at the faint end. There- 
fore, although different fitting techniques are at play, 
we have probed th e CO emission to the same depth as 
IBlanc et al.l (|2009[ ) , but at a much higher confidence level 
(2cr VS. la). 

We now investigate the impact of imposing threshold 
values to our gas maps. As can be seen in FiguresHHH the 
3-cr threshold has negligible impact on the SFR maps. In 
fact, varying the SFR threshold from 2a to 4ct changes 
the best-fit slopes minimally (within 0.03). This remains 
true at all scales. Most of the data censoring occurs on 
the gas maps. We test the reliability of our fits by varying 
the CO map threshold from la to 5a, and the results 
are tabulated in Table [3] and Table ID As can be seen 
from the various columns, decreasing the threshold of the 
CO map fiattens the slope 7H2 at all scales, suggesting 
that the inclusion of noisy data has adverse effect on the 
determination of the S-K law. With future deeper CO 
maps, our current 1-a limits may become better detected, 
and enable a more accurate determination of the slope of 
the S-K Law at low CO masses. Using our current well 
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detected data, with thresholds 3-<7 or higher, the slope 
of the molecular S-K law is clearly super-linear. 

Finally, one should be aware of the fact that because 
brightnesses at different wavelengths will often tend to 
correlate, causal relationships may be masked. In this 
work, this richness effect has been alleviated significantly 
when the underlying diffuse emission from dust and stars 
has been removed, and is thus of much less importance 
than in the study of BOS (one may notice from Fig- 
ure [3] that the background-subtracted images no longer 
appear disk-like). The absence of a clear correlation be- 
tween SFR and atomic hydrogen (especially when the 
local background component is subtracted, see Figure HI- 
d) is an implication that our analyses are not dominated 
by this effect. 

Further progress on this problem of how to measure 
SFRs locally within galaxies will require a careful mod- 
eling of the diffusion timescales of stellar populations in 
galaxies as they age, and the analysis of both the direct 
stellar light (UV, Ha) and the dust-reprocessed light (IR) 
from these evolving populations. 

6. SUMMARY 

We have presented a case study of the sub-kpc 
Schmidt-Kennicutt law of star formation in two nearby 
galaxies, M51a (NGC 5194) and NGC 3521, using our 
new CARMA + NR045 combined CO (J=l-0) maps of 
both galaxies with archival SINGS Ha, 24 fim, THINGS 
H I and GALEX FUV data. Employing a new approach 
of local background dete r minat ion, w e use the same data 
sets as IKennicutt et all (|2007D and IBigiel et al.l (|2008[ ) 
and quantitatively show that the contrasting results of 
the two studies on whether the molecular S-K law is 
super-linear or linear is a result of removing or preserving 
the local background, this being the sum of the diffuse 
emission in galaxy disks and bulges from stars and dust. 
We argue that in order to plausibly derive SFR maps of 
nearby galaxies, the local mid-IR and FUV background 
should be subtracted carefully, in order to isolate the 
currently star-forming regions. 
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Fig. 1.— The KPNO 2.1-m Ha, Spitzer MIPS 24 /<m, CARMA CO (J=l-0) and GALEX FUV images of NGC 3521. North is up, East 
is left. The scale of each panel is 4.0' x 4.6', corresponding to 9.3x10.7 kpc^. 
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Fig. 2. — Comparison of calculated molecular hydrogen surface densities in M51a between our CARMA + NR045 and BIMA SONG CO 
measurements at the 230 pc (5.8") resolution. The 3-a sensitivity limits of the SONG and our data at this resolution (18 and 9.9 Mq pc~^, 
respectively) are indicated by the dotted lines. The solid line indicates the expected unit slope in the case of perfect consistency. It is seen 
that the new CARMA + NR045 data recover more CO flux at the faint end, and is deeper by a factor of ~2 than the BIMA SONG data. 
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Fig. 3. — The products of our local background removal strategy applied on the 24 images of M51a (top row) and NGC 3521 
(bottom row). In each row, the compact emission component of the original image (first column) is identified and removed by adopting 
the algorithm of Hllphot, then the resultant image (second column) is smoothed by replacing each pixel with the median of a moving box 
so that no sharp features exist. The background image is thus created (third column) and subtracted from the original image, so that the 
final background-free image (last column) is prepared for our analysis. All images for each galaxy share common logarithmetic grey scales 
displaying surface brightness in units of MJy st~^. 
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Fig. 4. — The correlations for M51a between derived piiysical quantities at 750 pc resolution (at a distance of 8.00 Mpc). For the left 
panel in each pair the local backgrounds in the Ho, 24 jum and FUV images are not removed (denoted by "BG+"), but are removed in 
the right panel ("BG— "). The fitted slopes are indicated at the bottom-right of each panel, [a] the molecular-only S-K law with SFR 
derived from Ha-|-24 /im; [b] the molecular-only S-K law with SFR derived from FUV-l-24 /im; [c] the correlation between SFRs derived 
from Hq+24 fim and FUV-f 24 ^m, respectively; [d] the relation of Ha+24 fini SFR vs. H I surface density; [e] the total hydrogen S-K law 
with SFR derived from Ha-)-24 ^m; [f] the total hydrogen S-K law with SFR derived from FUV-l-24 fim. In each panel, the horizontal and 
vertical dotted lines locate the 3-o" sensitivity limit of the two involved quantities. Note that in general, the sensitivity limit of Sgpp)^ is 
not the limiting factor for our fitted range, while Sgas is. Those data points with sufficient S/N ratio for linear bisector fittings are shown 
as solid bla<;k dots, in contrast to the low S/N light grey dots. In panel [e], the solid line corresponds to the fit using points where CO, 
H I and SFR are independently above 3a (best-fit power index 7h = 1.63 ± 0.08), while the dashed line results from requiring the total 
hydrogen surfa<;e density and SFR to be above 3a (best-fit 7h = 1-42 ± 0.06). 
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Fig. 5. — The scaling relations of star formation for M51a at the highest available resolution (230 pc). The local backgrounds in Ha, 
24 fim and FUV images have been removed, following our strategy. The symbols, notes and involved physical quantities are identical to 
Figure [4] 
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Fig. 6. — The scaling relations of star formation for NGC 3521 at the highest available resolution for H I data (320 pc). The local 
backgrounds in Ha, 24 /im and FUV images have been removed, following our strategy. The symbols, notes and involved physical 
quantities are identical to Figure |4] All quantities are compared inside the FoV of CARMA CO data, which is ~ 30% of the optical area 
defined by D25 used by BOS. 
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Fig. 7. — Best-fit parameters as a function of apparent resolution (5): [a] tfie power index of molecular-only S-K law; [b] tiie r.m.s. 
dispersion of the data about the best-fit molecular S-K law in unit of dcx; [c] the power index of the total-hydrogen S-K law (requiring 3cr 
or better detection in both H i and CO maps); [d] the power index of FUV-l-24 /im SFR vs. Ha-|-24 fj,m SFR correlation. Note that from 
[a] to [c] the SFR is derived from Ha-t-24 //m luminosity. 
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Fig. 8. — The power index (7H2) and r.m.s. dispersion ((TH2) of the best-fit molecular S-K law in M51a and NGC 3521, plotted as a 
function of de-projected resolutions (S^p). Both quantities match remarkably well at a given S^p, suggesting a shared behavior of sub-kpc 
molecular S-K law in them. Fitting to the data points of M51a only, we find the logarithmic function = —1.09 log [5dp/kpc] -|- 1.36 
and the linear relation (TH2 = —0.19 [<5(jp/kpc] + 0.67 are able to depict the behavior of these two quantities in both galaxies remarkably 
well, as over-plotted with solid lines in [a] and [b], respectively. 
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Fig. 9. — The compact-to-diffuse ratio of emission in: [a] Ha, [b] 24 fj,m and [c] FUV, ploted as a function of Sh2- Only >3<7 pixels are 
taken into account. For both M51a and NGC 3521, the projected resolution examined is 400 pc. Note the broader range of variation in 24 
/im than in Ha and FUV. 
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TABLE 1 

Best-fit parameters of the sub-kpc molecular S-K law in M51a (S-cr thresholds are applied), expressed as SsFR(Ha+24pm) = ^ ^Yi2 ^'^^ 

an r.m.s. dispersion <th2 about the best-fit power scaling. 



Resolution (kpc) 


Th, 


A 




0.25 


1.86±0.03 


-5.49±0.05 


0.60 


0.30 


1.82±0.03 


-5.36±0.06 


0.60 


0.40 


1.65±0.04 


-4.97±0.06 


0.56 


0.50 


1.52±0.04 


-4.71±0.07 


0.55 


0.60 


1.49±0.05 


-4.62±0.07 


0.52 


0.70 


1.36±0.05 


-4.36±0.08 


0.50 


0.80 


1.32±0.05 


-4.24±0.07 


0.48 


0.90 


1.26±0.05 


-4.14±0.08 


0.44 


1.00 


1.21±0.06 


-4.08±0.09 


0.45 



TABLE 2 

Best-fit parameters of the sub-kpc molecular S-K law in NGC 3521. Symbols are the same as for M51a. 



Resolution (kpc) 


7H, 


A 




0.25 


1.41±0.06 


-4.48±0.08 


0.50 


0.30 


1.40±0.06 


-4.43±0.09 


0.45 


0.40 


1.35±0.09 


-4.32±0.12 


0.42 


0.50 


1.33±0.12 


-4.25±0.16 


0.39 


0.60 


1.19±0.11 


-4.05±0.15 


0.30 


0.70 


1.24±0.16 


-4.12±0.21 


0.27 



TABLE 3 

Best-fit slopes of the moleculaj: S-K law when different thresholds are applied but the SFR (traced by Ha-t-24 /Ltm) threshold is fixed 

at 3-cr: the case of M51a. 



Resolution (kpc) 


1-0- 


2-0- 


3-0- 


4-(j 


5-cr 


0.25 


1.60 


1.77 


1.86 


1.97 


2.02 


0.30 


1.51 


1.67 


1.82 


1.88 


1.97 


0.40 


1.41 


1.54 


1.65 


1.72 


1.75 


0.50 


1.33 


1.44 


1.52 


1.60 


1.66 


0.60 


1.32 


1.43 


1.49 


1.56 


1.63 


0.70 


1.20 


1.27 


1.36 


1.46 


1.48 


0.80 


1.21 


1.28 


1.32 


1.38 


1.48 


0.90 


1.10 


1.18 


1.26 


1.29 


1.32 


1.00 


1.06 


1.15 


1.21 


1.30 


1.33 



TABLE 4 

Best-fit slopes of the molecular S-K law when different Shj thresholds are applied but the SFR (traced by Ha-f 24 //m) threshold is fixed 

at 3-0-: the case of NGC 3521. 



Resolution (kpc) 1-<t 


2-<T 


3-0- 


4-0- 


5-0- 


0.25 


1.20 


1.31 


1.41 


1.47 


1.54 


0.30 


1.19 


1.34 


1.40 


1.44 


1.49 


0.40 


1.23 


1.26 


1.35 


1.35 


1.37 


0.50 


1.22 


1.22 


1.33 


1.33 


1.33 


0.60 


1.19 


1.19 


1.19 


1.26 


1.26 


0.70 


1.24 


1.24 


1.24 


1.24 


1.35 



